Theory of Flux-Flow Resistivity near Hc2 for s-wave Type-II Superconductors 
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This paper presents a microscopic calculation of the flux- flow resistivity p/ for s-wave type-II 
superconductors with arbitrary impurity concentrations near the upper critical field Hc2- It is 
found that, as the mean free path I becomes longer, p/ increases gradually from the dirty-limit 
result of Thompson [Phys. Rev. Bl, 327 (1970)] and Takayama and Ebisawa [Prog. Theor. Phys. 
44, 1450 (1970)]. The limiting behaviors suggest that pf{H) at low temperatures may change from 
convex downward to upward as I increases, thus deviating substantially from the linear dependence 
Pf(xH/Hc2 predicted by the Bardeen-Stephen theory [Phys. Rev. 140, A1197 (1965)]. 



PACS numbers: 74.25.0p, 74.25.Fy, 72.10.] 
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Kim et al. attributed finite resistivity observed 
in type-II superconductors to the motion of flux lines, 
calling it "flux-flow resistivity" p/. Whereas the idea 
has been accepted widely, we still have poor quantita- 
tive understanding of the phenomenon. The early phe- 
nomenological theories based on single-flux considera- 
tions d, 1^ Q cannot explain the steep decrease of p/ 
observed near Hc2 Q- This region, where microscopic 
calculations may be performed most easily, has been a 
subject of later theoretical works [aSSSl- We thereby 
have a quantitative theory on the flux-flow resistivity 
near Hc2 0) HI j but its validity is still restricted to the 
dirty limit. Reasonable agreements with the theory have 
been reported later in a couple of experiments on mi- 
crowave surface resistivity ;9i and flux-flow resistivity 
^d\. However, the latter experiment found a small dis- 
crepancy that p/ near Tc is larger than the theoretical 
prediction 0jjai which was attributed by Larkin and 
Ovchinnikov to the extra phonon scattering effective 
at finite temperatures. 

With these backgrounds, this paper provides a quan- 
titative theory on the flux-flow resistivity near Hc2 ap- 
plicable to arbitrary impurity concentrations. I thereby 
hope to establish the domain in which the dirty-limit the- 
ory BBl is valid, and look into whether the experiments 
[^I1(T|] may be explained by impurity scattering alone. Al- 
though a similar study was carried out by Ovchinnikov 
|l2l [. the behaviors of p/ at intermediate impurity con- 
centrations have not been clarified explicitly. 

I consider the s-wave pairing with an isotropic Fermi 
surface and the s-wave impurity scattering in an external 
magnetic field H || z. I calculate the complex conductivity 
cf{uj) below microwave frequencies uj using the quasiclas- 
sical equations of superconductivityllj] , and finally put 
UJ ^ 0. I use the notation of Ref. [ij but leave the Hall 
terms out of consideration. Hence I start from the same 
equations as Eschrig et al. in clarifying the motion 
of a single flux line within the linear-response regime. 

The vector potential without perturbation is given by 
A(r) = Bxy+ A(r), where B is the average flux density 
and A expresses the spatially varying part of the mag- 



netic field satisfying J V x A dr = 0. The corresponding 
retarded quasiclassical Green's functions and can 
be obtained as Ref. with the replacement of the Mat- 
subara frequency e„ by —ie. I adopt the units where the 
energy, length, magnetic field, and electric field are mea- 
sured by the zero-temperature energy gap A(0) at H = 0, 
the coherence length = ^''^f /^(O) with vp the Fermi ve- 
locity, i?o = 0o/27r^Q with (j)o = hc/2e the flux quantum, 
and EQ = hv-p I eS^Q, respectively. I also put h — k^ — \. 

Now, consider the response to a spatially uniform but 
time-dependent perturbation JAe^*'^* = 5Ee~''^*/*^ with 
(5E_LH. A straightforward calculation based on Eq. (71) 
of Ref. Q leads to the following equation for the flrst- 
order response 5f^ — 6f^{e,'k,L0,r): 
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Here t is the relaxation time in the Born approxima- 
tion, (• ■ ■ ) denotes the Fermi-surface average satisfying 
(1) — 1, and A(r) and <5A(u;,r) are the pair potential 
and its first-order response, respectively. The unit vector 
V = k specifies a point on the spherical Fermi surface, 
.f±=/^(£±,k,r) and 5± =3^(£±, k, r) with e±=£±w/2, 
and the superscript ''' denotes simultaneous operations of 
complex conjugation and (e, k, lj)^ (— e, — k, — w), e.g., 
8(^'^{e, k, g;, r) = 8(^*{—e, — k, —lo, r). Finally, the nor- 
malization condition [ij] enables us to write in terms 
of 8]^ as 

'53^ = - ( <5./«^ + f-' Sn/{ gl + g^). (lb) 

Equation ^ determines the retarded functions Sf^ and 
Sg^ for given SE and <5A. In addition, the advanced func- 
tions are obtained directly by using Eq. (72) of Ref. |3 
as Sf^{e, k, uj, r) — Sf^{~e, — k, uj, r) and Sg^{e, k, uj, r) = 
~Sg^*is,k,-iu,r). 
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As for the Keldysh functions, I write them following 
Eschrig et al. as 5g^ = Sg^cf)- - (f>+Sg^ + Sg"" and 
5f^ = 6f^(j)_-(j)+5f^+6f'' with 0± = tanh(e±/2r). Then 
a simplified equation results for (5g* = (5.g^(£, k, a;, r) as 
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It also follows from the normalization condition that 
df^ = Sf^{e, k, uj, r) is given in terms of 6g^ by 

Sr = -{ Sg^^ + 5g- ) / ( 5? - .9^ ) • (2b) 
Finally, the pair potential (5A is determined by 



6A 



de{5f''cj,^~5f^cj,+ + 6n. (3) 



where -A''(0) is the density of states per spin at the Fermi 
level, Vq > is the s-wave pairing interaction, and Sc 
is the cut-off energy. Once the solution to Eqs. Q-® 
is obtained self-consistently, the transport current 5] — 
(5j(a;, r) is calculated by 
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V+ + -55")). (4) 



The corrections caused by 5] to the magnetic field and 
the charge density may be neglected safely for the rele- 
vant weak-coupling superconductors. Equations 
are still exact within the linear-response regime and form 
a basic starting point to calculate complex conductivity 
of type-II superconductors at arbitrary magnetic fields. 

I now concentrate on the region near Hc2 and solve 
Eqs. Q-® by expanding every quantity up to second or- 
der in A(r) as /^^/^d), gi^^l-i/RtD/^CDt, A = k^^\ 

with (^R^") =0 as seen from Eq. Ijlb|l . The superscript 
in 5f and (5A will be dropped as it causes no confusions. 

The zeroth-order quantity Sg'^^^^ is obtained easily 
from Eq. H2a|l as 
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When put into Eq. Q , this expression yields the normal- 
state Drude conductivity (T„, as it should. 

Next, Eq. (|la|l is simplified for the first-order 5f^ into 
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where a^^ and a are creation and annihilation operators 
satisfying [a, at] = 1 m 113, and {9, if) are the polar 
angles of v. Equation © can be solved with the Landau- 
level expansion (LLX) method [T^ by expanding SA and 
(5/^ in periodic basis functions of the vortex lattice as 

oo 

SAiiu, r) = VvJ2 S^N{i^) ^AjVq(r) , (7a) 

oo oo 
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where N denotes the Landau level, q is an arbitrary 
chosen magnetic Bloch vector, V is the volume of the 
system, and ■(/'Arq(r) satisfies ai^Nq = VN'^N-iq and 
0''^'>pNq = V-^+lV'-ZV+iq- The quantities without pertur- 
bations are expanded similarly. Near Hc2, A(r) may well 
be approximated using only the = level as A(r) = 
VV Aq V'oq(r). The coefficients Aq and /,^^ have already 
been obtained in Ref. 0, satisfying f^j^ — SmNAof^ 
with f^*{e,e) = f^{-e,0). It then follows from Eqs. @ 
and (j^l that Sf^^jy and JAjv may be written as 

6An = Ao(5a,i5Ai , (8a) 



= Ao (^Sm.N-l 5/^+ + S^.N+l <5/^_) . (8b) 
Equation is thereby transformed into (/i = ±) 
E MNN'6f^,^ = (/^+ + f^-)SE^ sin e/Au; 

N' 

+6niS^+ (<5Ai + (^/T+)/2r) , (9) 

where E± = E^±iEy, /^± = /||(e±), and 

Mnn' = -HSnn' + (3VN + 15n^n'-i - PVn Sn,n'+i , 

(10a) 

with e = £+i/2T and /5ee sin 61/2^2. To solve Eq. ©, 
let us write the inverse of the matrix A4 as 



A^' = {M-^)nn' 
Then Sf^^ is obtained formally as 



(10b) 



6f^^ = K^SE^/u; + Sf,+Kl, (^SAi + {6jf^)/2T) , (10c) 



where K?; is defined by 
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with /jy already obtained as 01 

/^ = n/(l-«)/2r). (lOe) 



Taking the angle average of Eq. (|10c(l , we obtain 
(iff) 6E+ , {KD 



1+/ 



1 - {KI)/2t uj 1 - {KI)/2t 
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Equation pU|l determines Sf^^ efficiently and fixes the 
retarded response Sf^. 

We next consider Sf^ and substitute Eqs. l|7b(l . and 
(pbjl into Eq. l|2b)l with = -g^^l. Expanding (5/*" as 
Eq. Ij7b|l . we find that the coefficient (5/^^ can also be 
written as Eq. I|8b|l . The corresponding Sff^^ {fi = ±) is 
obtained easily as 



with {fi = ±) 
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The quantity = qj._j_i n can be obtained simi- 
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larly from Eq. H2a|) as 

^(e,c.) = Agr^,(e,c.)^M^^<5£;^ , 
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Now we are ready to calculate the pair potential SA. 
Let us substitute the above results for Sf^ and Sf^ into 
Eq. 10). Then a straightforward calculation yields the 
expression for SAi defined by Eqs. H7a|l and (|8a|l as 



SAi{lu) = D{uj) 
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In deriving the formula, use has been made of the sym- 
metrics: K};{e) = Klf{-e), Kf}{s,uj) = K^ie,-Lu) = 
K^*(-e,Lu), and KUe^co) = (-D'^'^Kl 
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de 
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Kl{e) — > i/e as e 



I also substituted the result l/7V(0)Vb 



at T = = and put Ec ~^ oo, noting 
oo, etc. One can see easily that 



Eq. lfT2H|) satisfies £>(-w) = D*(w). 

The second-order quantities (Sg^'^'^' can be calculated 
similarly by expanding them as 

5<7^-(^) (e, k, r) = ^ ^ 8gl^'^ (e, ^, c.) e""-+^^- , 



rn K 



(13) 

where K is a reciprocal lattice vector of the magnetic 
Brillouin zone 17]. Since (5E is spatially uniform, we only 
need the K = component within the linear-response 
regime. It also follows from Eq. Q that m ^ ±1 are 
irrelevant for the current density. Let us substitute the 
expansions for and 5f^ into Eq. I|lbp with g^ = 1, 
multiplies it by e^*'''/27ry, and perform integrations over 
(r, ip). We thereby obtain g^ 



-ym=±l,K=0 



Alrlis,.)'-^ 



(14a) 
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Substituting Eqs. (fT!^ - lfT5|) into Eq. ^ and noting 
r^{s,u) = r«*(-e,-^) and T-{e,Lo) - Tt^{~e,Lo) = 

r^j*(— e, ~ll>), we finally obtain Sj^"^ = ji^''± ijy'^'' as 



<5ji'^(c.)=af)(a;)5i?±, 



(16a) 



with 



af\oj)^-Al 



3N{0)v F 
2(7 



;E([r>,^)-r;i(£,c.)]sin(?). (16b) 



("2") 

Thus (Ty is isotropic, as expected for any second-order 

corrections, and satisfies {—uj)^(j^j^''* (uj). 

Equation (|16b|) forms the main result of the paper, 
which is not only exact but also convenient for nu- 
merical calculations. The key quantities are 's de- 
fined by Eq. Ijl0b|) with Eq. (|10a|) . as may be seen from 
Eqs. iTOd|) .^ fTUe^ . ((TTb|) .^ (pb|) . fMbj) . and (pb^ . sat- 
isfying = (^_]^)JV+Af x^,. An efficient algorithm 
to calculate them is obtained as follows. Let us de- 
fine Pjv (T^n) for = 0,1,2,- •• as the determinant 
of the submatrix obtained by removing (retaining) the 
first N rows and columns of the tridiagonal matrix M , 
Eq. (|lUa|l . They satisfy Pjv-i = -HVn + P^NVn+i 
and T>N+i — —iS'Dn+P'^N'Dm-i with Vi = — ie and 
2?o = 1, as can be shown with standard matrix ma- 
nipulations. Then for A^' < is obtained, by 
also using standard techniques to invert a matrix, as 
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FIG. 1: The slope a of the flux-flow resistivity at Hc2 as a 
function of the reduced temperature T /Tc for several values 
of 5e/^ = 1 with KGL = 50. The curve TTE denotes the 

prediction of the dirty-limit theory 0, • 



This algorithm can be put into a more convenient form 
in terms of TZn = T^n+i/T^n and TZn = T^n-i/T^n as 
follows. Let us calculate TZn and 7?.jv by 



Un-i = i-ie + P^NTlN) 
Un+i ^ i-ie + P^nUn) 



-1 

-1 



^1 = 



i/e, 



substantially over the prediction of the dirty-limit the- 
ory. This result near Hc2 also suggests that pf{H) at 
low temperatures may change from convex downward to 
upward as I increases and may not be fit quantitatively 
by the Bardeen-Stephen theory p/ oc H/Hc2 0- The 
slope has also been found to become steeper as kgl ap- 
proaches l/-\/2, due to the increase of the coefficient of 
Al(xHc2-H, reaching a = 9.57 (2.59) at T = Tc (r = 0) 
for KGL — 1 and S,e/1 = 50. This fact indicates the neces- 
sity of correctly identifying the material parameters, such 
as kgl, and I, for any detailed comparisons between 
the theory and experiments on the flux-flow resistivity. 
Besides, a careful experiment will be required, especially 
near to determine the slope a which may change ap- 
preciably near Hc2 0- 

In summary, this paper has developed a reliable and 
efficient method to calculate the flux-flow resistivity near 
Hc2 over all impurity concentrations and clarified large 
dependence of p/ on both I and kql- 

This research is supported by Grant-in- Aid for Scien- 
tific Research from the Ministry of Education, Culture, 
Sports, Science, and Technology of Japan. 
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for an appropriately chosen large iVcut- 
N'<N is obtained by 



Then K^' for 
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One can check the convergence by increasing iVcut- It 
turns out that Ncut = 2 is sufRcient both near Tc and in 
the dirty limit where an analytic calculation is also possi- 
ble. One can thereby reproduce the formula by Thomp- 
son and Takayama and Ebisawa ^ which satisfies 

(2) 

— > 0) = 0. In contrast, A^cut ^ 1000 is required 
in the clean limit at low temperatures. 

Using Eq. Ijl6b|l and taking the limit w ^ 0, I have 
calculated the initial slope of the flux-flow resistivity: 



H dpf 
Pn dH 



(18) 



for various impurity concentrations and various values 
of the Ginzburg-Landau parameter kql- The main H 
dependence in Eq. I)16b(l lies in AqOc Hq2 — H , which can 
be obtained accurately following Ref. 16. 

Figure 1 displays the slope a for several values of 
= 1/2ttTcT with KQL = 50. The curve TTE de- 
notes the prediction of the dirty-limit theory by Thomp- 
son and Takayama and Ebisawa Marked mean- 
free-path dependence is clearly seen. In fact, the slope 
at r = Tc (r = 0) decreases from 4.99 (1.72) in the dirty 
limit to 0.89 (0.73) at = 4.0. Thus, nonlocal ef- 
fects in clean systems tend to increase the resistivity 
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